Spintronic Properties of Zigzag-Edged Triangular Graphene Flakes 
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We investigate quantum transport properties of triangular graphene flakes with zigzag edges 
by using first principles calculations. Triangular graphene flakes have large magnetic moments 
which vary with the number of hydrogen atoms terminating its edge atoms and scale with its 
size. Electronic transmission and current-voltage characteristics of these flakes, when contacted 
with metallic electrodes, reveal spin valve and remarkable rectification features. The transition 
from ferromagnetic to antiferromagnetic state under bias voltage can, however, terminate the spin 
polarizing effects for specific flakes. Geometry and size dependent transport properties of graphene 
flakes may be crucial for spintronic nanodevice applications. 

PACS numbers: 73.63.-b, 72.25.-b, 75.75.+a 



I. INTRODUCTION 

Graphene and graphene based nanostructures are fo- 
cus of intensive research activities due to their impressive 
material properties [THE] and promising application po- 
tential in novel electronic devices. In particular 
edge-localized spin polarizations found in graphene rib- 
bons [E2QJ], flakes [TSUI?] , and at defect sites PUCES] in- 
troduce magnetic properties that can be utilized for spin- 
tronic applications. Recent studies have also revealed the 
ferromagnetic ground state of graphene nanodots|20j. tri- 
angular shaped graphene fragments [2"Tr[2~4] and graphene 
domains on 2D hydrocarbons [25] and the possibility of 
observing spin polarized current voltage characteristics 
of such graphene flakes. 

While pristine graphene provides high carrier mobility 
and ambipolar behavior, semiconductor nanoscale ma- 
terials having tunable bandgap are more desirable from 
the perspective of potential nanoelectronics applications. 
In this context, recent efforts have been devoted to pre- 
cise controlling electronic and magnetic properties of 
graphene sheets by functionalization via adatom adsorp- 
tion. The synthesis of a 2D hydrocarbon in honeycomb 
structure, namely graphane . [26H29] is one of the succes- 
full example for such functionalization. Very recently, 
we reported the possibility of obtaining tunable bandgap 
and magnetization through dehydrogenation of domains 
on 2D graphane and graphane nanoribbons.|25| [30] Sta- 
bility and electronic properties of graphene flakes uni- 



formly functionalized by methyl (CH 3 ), phenyl (QqR 5 ) 
and nitrophenyl (C6H4NO2) groups were also discussed 
earlier. [3IH33] 

Recent experimental observations [34-38 and theoret- 
ical studies [3"§H4T] show that the electronic, magnetic 
and conductance properties of graphitic fragments can 
be changed significantly upon the termination of their 
edges. In addition to purely zigzag and armchair edged 
graphene, experimental verification of the existence of 
alternating series of zigzag and armchair segments at 
the edges and energetics of reconstructions have also 
been reported. |42 Originating from the antiferromag- 
netic ground state of zigzag edges, adatom[15] and 
topology[32] dependent trends in electronic properties of 
rectangular flakes have also investigated. 

In this work, we study graphene flakes having equi- 
lateral triangular shapes with zigzag edges (n-TGF), 
where n denotes the number of edge hexagonal cells in 
one side of the triangle. The flakes have been consid- 
ered as bare (C„2 +4n+1 ), each edge atom being satu- 
rated with one (C n 2 +4n+1 H3„ + 3) or two hydrogen atoms 
(C n 2 +4n+1 H6n+6)- We find that these flakes have large 
spin magnetic moment values of 4(n — 1), (n — 1), and 
2(ra — 1), respectively, in units of \ib- When these TGFs 
have been contacted with thin metallic electrodes we cal- 
culate that the current running through them gets both 
spin-polarized and rectified. 



II. CALCULATION METHODS 
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Optimization of geometrical structures of triangular 
graphene flakes and calculations of their magnetic and 
electric properties are performed by using the software 
package Atomistix ToolKit (ATK) [H] based on density 
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FIG. 1: (Color online) Atomic structure and isosurface of 
charge density difference of spin-up (t) and spin-down (J,) 
states for 4- (upper row) and 5-triangle (lower row) graphene 
flakes: Bare, singly- and doubly-hydrogenated edges. Calcu- 
lated net magnetic moments of the flakes are given in terms 
of Bohr magneton (/xb)- Difference charge density of spin- up 
and spin-down states is shown by red (dark) and blue (light) 
isosurface, respectively. 




FIG. 2: (Color online) Electrode-device geometry and con- 
vention for forward and reverse bias applied to triangular 
graphene flakes (TGFs). 



the TGFs is not fixed and energy is minimized with re- 
spect to the electrochemical potentials of the electrodes 
at each voltage increment. In order to achieve conver- 
gence of the electronic states with increasing voltages 
to the desired level of accuracy, calculations performed 
within the bias window (from -1 to +1 V) in steps of 0.005 
V. By using the carbon chains attached to the TGF, it 
is ensured that the screening occurs in the device region. 



functional theory (DFT). The spin-dependent exchange- 
correlation potential is approximated within the general- 
ized gradient approximation. The criteria of convergence 
used for total energy and Hcllman-Feynman forces were 
10~ 4 eV and 0.005 eV/A, respectively. The electrostatic 
potentials were determined on a real-space grid with a 
mesh cutoff energy of 150 Ry. Double-zeta-polarized ba- 
sis sets of local numerical orbitals were employed to in- 
crease the accuracy of our calculations. 

For determination of quantum transport properties of 
the electrode-TGF-electrode system, ATK use nonequi- 
librium Green's function formalism. Transport calcu- 
lations are performed with the Brillouin zone sampled 
with (1,1,51) points within the Monkhorst-Pack k-point 
sampling scheme. The current through the TGFs is de- 
termined by summing the transmission probabilities for 
electron states from one electrode to another within the 
energy window [i^ — \ir = |eV|, where /i^ is the 

electrochemical potential of the left (right) electrode un- 
der the applied bias V . Therefore the spin dependent 
current is given by the formula 

I*(V) = G r R T a (E,V)dE (1) 

J f-L 

where Go = e 2 Jh is the quantum conductance unit and 
T a (E, V) is the quantum mechanical transmission proba- 
blity for electrons with spin state a and energy E. During 
the self-consistent calculation of I-V spectrum charge on 



III. RESULTS 

A. Atomic and Magnetic Ground State Properties 
of TGFs 

First, we have performed geometry optimizations of 
TGFs and determined the spin polarized charge den- 
sity of the optimized structures. Due to the particu- 
lar shape of the TGFs spin-relaxed calculation leads to 
lower ground state energies compared to those of spin- 
unrelaxed calculation; here we have found that ferromag- 
netically ordered spin accumulation at the edges gives a 
nonzero magnetic moment to the flake. Apart from minor 
bond contractions of the edge atoms TGFs preserve their 
regular hexagonal structure even when the flakes are not 
hydrogenated. Saturating the edge carbon atoms with 
either one or two hydrogen atoms considerably modifies 
the magnetic ground state, by altering the total magnetic 
moment. 

Graphene, with its hexagonal lattice structure result- 
ing from sp 2 -type hybridization of carbon atoms, is a 
planar ir conjugated system. It can be viewed as made 
up of A- and B-sublattices of carbon atoms. Repulsive 
case of the Lieb's theorem [35] reveals uniqueness of the 
ground state and provides a simple formula for calculat- 
ing the magnetic ground state of such bipartite systems. 
According to the rule provided by the theorem, total net 
spin magnetization of a graphene structure is given by 
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FIG. 3: (Color online) (a) Energy level spectra within ±1 eV range of Fermi level (E F ), the HOMO-LUMO gap (A) and 
isosurface of HOMO and LUMO orbitals. (b-c) Calculated I-V curves for hydrogenated and bare triangular graphene flakes 
(TGFs). Results of 4-TGF and 5-TGF are presented, respectively. Spin-up (t) and spin-down (4-) currents are shown by red 
(dark) and blue (light) lines, respectively. Solid and dashed lines denote forward and reverse bias calculations, respectively. 
Transmission curves of spin-up and spin-down under zero bias are also shown by insets. Transmission spectra of singly-H and 
bare TGFs are plotted up to maximum value of 0.03 and 0.4, respectively. Fermi levels are set to zero. 



M= -\N A - N B \gp B (2) 

where g ~ 2 for electron, /is is Bohr magneton, and 
Na and Nb denote the number of carbon atoms in A- and 
B-sublattice, respectively. Each carbon atom in graphene 
is connected to the nearest neighbors by three covalent 
bonds, while leaving behind a j^-orbital electron con- 
tributing to the spin magnetic moment. Such electrons 
in each sublattice have opposite spin states and make 
spin paired 7r-bonds, thus total spin magnetic moment 
of the system is zero unless a difference is created in the 
numbers of atoms of the A- and B-sublattices. In some 
graphene flakes, such as those having equilateral trian- 
gular shapes and zigzag edges, Na and Nb are different 
and leads to finite net spin magnetic moments (32 H3] • 
Moreover, in bare (not hydrogenated) flakes the edge car- 
bon atoms that make only two covalent bonds contribute 
to the spin moment with two nonbonding electrons. 

In Fig[l] we present the optimized structures and spin 
charge density difference (Ap = — p^) isosurface for 
bare, singly- and doubly-hydrogenated cases of n-TGFs 
(n=4, 5). In general, we see that each carbon atom has 



an induced spin imbalance opposite to its nearest neigh- 
bors, however the spin polarization of hydrogen atoms is 
negligible. 

Calculated total spin magnetic moments of the struc- 
tures using DFT are in integer multiples of ps and ei- 
ther verifies the results of Lieb's theorem or can be un- 
derstood by simple modifications of it. In the flakes 
with singly-hydrogenated edges, all the carbon atoms 
are coordinated as reminiscent of infinite graphene, and 
pj Pb = N A — Nb = n — 1 in accordance with Lieb's the- 
orem. When the hydrogen atoms are removed from the 
flake (bare flake case) the 3n atoms of sublattice A and 
3 atoms of sublattice B each have an extra nonbonding 
electron that contributes to the magnetic moment, con- 
sequently giving p/ Pb = n— I + 3n — 3 = 4(n — 1). On 
the other hand, in doubly-hydrogenated flakes, all the 
edge carbon atoms have sp 3 hybridized electrons with 
vanishing contribution to spin magnetic moment. Thus, 
in this case p/pb = \ n — 1 — 3n + 3| = 2(n — 1), where 
the direction of magnetic moment is reversed, and the 
net moment of the system increases with respect to the 
singly-hydrogenated case. Tunability of the spin mag- 
netic moments of triangular graphene flakes through de- 
gree of hydrogenation is an interesting feature that may 
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FIG. 4: (Color online) Bias dependent transmission spectra of singly hydrogen passivated (singly-H) and bare TGFs for forward 
bias. Spin-up and spin-down transmissions are shown by red (dark) and blue (light) lines. Ef is set to zero. 



be utilized for nanodevice applications. Saturation of 
TGF edges by other atoms may give rise to similar mod- 
ifications in their electronic and magnetic structures. 



B. Transport Properties of TGFs 

Next we consider electrical conductance and I-V char- 
acteristics of TGFs when contacted with metallic elec- 
trodes. We use linear carbon chains as model elec- 
trodes. Carbon atomic chains which are known to be 
metallicfTBl |4"o] are expected to make reasonably good 
contacts with the flakes. As an alternative to carbon- 
chain electrodes, earlier the robustness of our conduc- 
tance calculations was also tested by using semi-infinite 
gold bar electrodes and consistent results were obtained 
for rectangular flakes. [15] At the contact sites the intro- 
duced hydrogen atoms are removed and connection be- 
tween electrode and flake carbon atoms is established 
through double bonds. Due to asymmetric form of 
device-electrode geometry one expects different current 
profiles for forward and reverse bias voltages. The con- 
vention we have adopted is schematically shown in Fig(2] 
where forward bias corresponds to flow of electrons from 
the left to the right electrode, i.e. the current is flowing 
from right to left electrode. 

For the ballistic current of electrons, the spin de- 
pendent transmission probability from one electrode to 
other strongly depends on the eigenstates of the TGF 
molecule. The molecular orbitals perfectly coupled to 
the electrodes behave as an open channel and provides 
propagation with minimal scattering through the TGF 
and hence the magnitude of the transmission coefficient 
is determined by electrode-device interaction. Eigen- 



states within the Ep±l eV energy window, indicat- 
ing the highest-occupied-molecular-orbital (HOMO) and 
lowest-unoccupied-molecular-orbital (LUMO) of singly- 
hydrogenated and bare n-TGFs (n=4,5) are presented 
in Fig{3] (a). It is obvious that the HOMO-LUMO gap 
(A) gets narrower upon the hydrogenation of TGFs. Re- 
sulting from the ferromagnetic ground state, within the 
Ep±l eV energy window up and down spin states are 
well-separated around the E^. In the energy window 
used for Fig{3] (a), all the filled levels are for up-spin 
states, whereas the unoccupied levels are for down spins. 
When the TGF molecule is connected to electrodes, how- 
ever, due to the chain-TGF interactions energy level spec- 
tra is changed and it is not easy to exactly distinguish 
the contribution of TGF and electrode states. 

In Fig|3](b-c), we show the I-V characteristics of 4- and 
5-TGFs for bare and singly-hydrogenated cases. Hydro- 
genated flakes have lower conductance and leads to at 
least an order of magnitude smaller currents, since the 
hydrogenation of the flake removes some of the states 
providing open channels in the flake. While the maxi- 
mum current in the calculated bias range is ~ 5.00 /iA 
for bare TGF, after hydrogenation it is reduced to ^0.35 
fj,A. Since electrode device coupling broadens the energy 
levels, and may shift them due to charging, even though 
there is no molecular state at Ep, the tails of HOMO and 
LUMO states can contribute to the transmission even at 
small voltages. Orbital character that changes upon the 
H termination of HOMO and LUMO states reveal the 
importance of the edge atoms in electron transport. In 
addition, there is a strong rectification of current for both 
singly-H and bare edge cases. In this electrode-device 
configuration forward current gets larger at a threshold 
bias of ^0.6 V. For small voltages, the characteristics 
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of the current flowing through the TGFs can be under- 
stood by zero bias transmissions that are given by insets 
in Fig|3] (b-c). In fact the self consistently calculated 
voltage dependent transmission spectra provides a bet- 
ter information regarding the I-V characteristics of the 
device. 

Depending on the applied voltage, triangular graphcne 
flakes display diverse properties. In the case of for- 
ward bias application to the singly-hydrogenated 4-TGF, 
the spin-up current is dominant up to 0.55 V, but af- 
ter this critical value, the spin-up and spin-down states 
are merged due to the transition from ferromagnetic to 
antifcrromagnetic state and thus spin polarizing prop- 
erty disappears. This behavior of I-V curves can be 
revealed from transmission spectrum. Even if the zero 
bias transmissions shown in insets in Figj3] (b-c) can 
explain I-V curves for small bias values, at finite bias 
voltages the transmission spectrum of the flake changes 
through broadening and/or shifts of the transmission 
peaks, which should be calculated self-consistently under 
nonequilibrium conditions. In Fig|4] we show the vari- 
ations of transmission spectrum for selected cases un- 
der incremental forward bias voltages. For the hydro- 
genated 4-TGF, while transmission of the up spin chan- 
nel decreases with increasing voltages, the transmission 
of down spins is increasing. Eventually both transmis- 
sion curves are merging at the vicinity of 0.55 V and 
hence the spin polarization of the current is ceasing. To- 
tal magnetic moment of the flake together with electrodes 
has a bias dependence, gradually decreasing from 2 /ig 
to zero at the merging point of the up and down spin cur- 
rents. In contrast, the down spin states of bare 4-TGF 



under forward bias get closer to the E p and become dom- 
inant carriers in the current. This explains how the spin 
polarization of the current is switched by applied bias. 
We can also establish a relation between I-V behavior 
and the corresponding transmission spectrum of bare 5- 
TGF. Since up and down spin transmission peaks do not 
show considerable variation within the energy window, 
the spin polarization of the currents is maintained up to 
1 eV. 



IV. CONCLUDING REMARKS 

In summary, we have investigated the electric, 
magnetic and transport properties of triangle shaped 
graphene flakes. We have found that in addition to 
the their ferromagnetic ground state, triangular graphene 
flakes show spin polarized and rectified current properties 
depending on edge saturation, flake size, bias voltage and 
bias direction. Diverse and spin dependent properties of 
graphene flakes depending on their shape, size and edge 
saturation keep the promise of variety of application in 
future nanospintronics. 
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